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Abstract. We study a general model of granular Brownian ratchet consisting of an 
asymmetric object moving on a line and surrounded by a two-dimensional granular 
gas, which in turn is coupled to an external random driving force. We discuss the two 
resulting Boltzmann equations describing the gas and the object in the dilute limit and 
obtain a closed system for the first few moments of the system velocity distributions. 
Predictions for the net ratchet drift, the variance of its velocity fluctuations and the 
transition rates in the Markovian limit, are compared to numerical simulations and a 
fair agreement is observed. 
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1. Introduction 

A Brownian ratchet is a system designed to extract work, usually in the form of a 
net drift or current, from a thermal bath. If the bath is at equilibrium, i.e. it is 
characterized by only one temperature, a ratcheting behavior is prevented by the second 
principle of thermodynamics. The system must be coupled to different baths at different 
temperatures, or to additional specific non-conservative (e.g. time-dependent) forces, 
in order to escape the consequences of the second principle. Moreover, to observe a net 
drift, spatial symmetry must also be broken [TJ [2j [3]. Dissipation of energy in inelastic 
collisions between macroscopic grains [I], breaks the time-reversal symmetry and leads 
to the introduction of suggestively simple models of inelastic Brownian ratchets, which 
are apparently coupled to only one thermal bath at a single temperature. It is known, 
for instance, that a granular object, surrounded by a stationary inelastic gas and 
characterized by a left-right asymmetry, presents a rectification of thermal fluctuations 
[5] [6] resulting in a net drift in a given direction: the asymmetry of the object can 
originate from its shape or its inelasticity profile (i.e. different inelasticities on different 
portions of the surface). Interestingly, the probability density function (pdf) P(V) of 
the object velocity results asymmetric when its mass M is of the same order or smaller 
than that of surrounding disks [TJ: such an effect is stronger the smaller the elasticity. In 
this study we intend to offer an analysis of a general model which includes all ingredients 
cited above, treating also the dynamics of the surrounding gas. This will make clear 
the conditions required to decouple the gas dynamics from that of the object, making 
the latter obey a closed Markovian master equation. In this limit we will obtain some 
general formula for the net drift of the object, its velocity variance and the transition 
rates of the Markov process, which in general do not satisfy detailed balance. These 
results compare very well with numerical simulations. The present article is organized as 
follows: in Secj2] we introduce the model and obtain the coupled equations of evolution 
for the probability distributions of the velocities of the ratchet and of the gas, in Secj3] 
instead of solving directly the former equations we consider the governing equations for 
the moments of the distributions, while in Sec. 0] we specialize our study to the case of 
an equilateral triangle, while in Sec [5] we conclude with a brief summary. Finally, we 
provide an appendix containing the necessary formulas for the coefficients entering the 
equations of Secj3] 

2. Theory 

Our 2D model consists of a rigid convex object of mass M, generally asymmetric, 
surrounded by a dilute gas of N hard disks of mass m and density n = N/ A where A is 
the area of the box. The surface of the object, of perimeter C, has a non- homogeneous 
inelasticity with a coefficient of restitution that depends on the point of contact on 
the surface. The object can only slide, without rotating, along the direction x. The 
collisions between two disks of the gas are dissipative with a coefficient of restitution add- 
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In this system the energy is not conserved and an external driving mechanism is needed 
to attain a stationary state. In order to maintain this steady state the gas is coupled to 
a thermal bath [8]: the gas particles between two collisions are subject to an external 
random force. The dynamics of the object is assumed not to couple directly with the 
thermostat, but only with the gas particles. For the sake of simplicity we assume that 
the gas is dilute and that Molecular Chaos is valid for object-disks collisions: this allow 
us to use the Direct Simulation Monte-Carlo (DSMC) algorithm to simulate the system 
dynamics [9]. 

After a binary collision, the velocities of the particles and the object, v and 

V respectively, can be obtained, from their pre-collisional values v 1 = (v' x , v' y ) and 

V = (V, 0), imposing the following conditions concerning a portion of the surface: 

MV + mv x = MV' + mv' x (1) 
u-jfe|| =v'-k\\ (2) 

{V - v) ■ k± = - a{9){V' - if) ■ k ± . (3) 

k\\ and k± are the unit vectors, parallel and perpendicular respectively, to the object 
surface in the collision point (see Fig. [T]). These can be expressed as fcii = (cos 9, sin 9) 
and k\\ = (sin 9, — cos 9) where 9 is the angle created by kn and the x axes, modulus 2n 
(Fig. CEJ). In this way the coefficient of restitution is a function of 9 and it is expressed as 
a(9). Equations ([1]) and ([2]) correspond to momentum conservation in the x direction 
and parallel to the surface, while the inelasticity takes part only in the reduction of 
the relative velocity expressed in Eq. (j3J). Solving these equations, we can write the 
post-collisional velocities as 

T w [I + a(9)]e 2 sm9 . tW . . „ , m 

V = V + i ^2 2n \i< ~ V) sm 6 ~ K cos 0} 
1 + e l sm 9 

, t 1 + sin6» , . , 

v x = v x - ! 2 \{v x -V)sm9- v CO8 0] 

1 + e z sm 9 

. [1 + a (9)] cos9 r . . . . . 

v » = v v + i -l 2 2 fl « -V')sm9- v' y cos 9} (4) 
1 + e z sm & 

where e 2 = m/M is the mass ratio. The collision between two disks labeled 1 and 2, 
instead, is given by 

_ 1 + add * * 
Ui = (ui2 • njn 

V2 = v' 2 + 1 (t7 12 . n)n, (5) 

with v\2 = Vi — V2 and n the unit vector in the direction joining the centers of the two 
disks. As external driving force acting on the disks, we choose the following heat bath: 
d ~ * 

77t— — = -mTv{t) + ^2mTT b ({t) (6) 

with ((t) a Gaussian white noise with (<f(£)) = and (Ci(t)Cj(t')) = — an d 

T = l/r b the drag coefficient. The model has been studied in 2D [8] and ID [10] and 
also without viscosity [11]. In the dilute gas limit, it is possible to describe the object 
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dynamics by means of a Boltzmann Equation (BE) for P(V, t) which can be written, if 
the object is convex, as [5] 

^ = J dV [W(V\V')P(V',t) - W(V'\V)P(V,t)] (7) 
where the rates for the collision object-disk is 

r2n /»+oo /•-f-oo 

W(V\V')= nCF(9)d9 / dv' x / dv' y <f>(v' x ,v' y ,t) 

Jo J— oo J — oo 

G[(V' - v 1 ) ■ k ± }{V' - v) ■ k±5tv - V 

with K>(9,a(9)) = (1 + e 2 sin 2 9) /[{I + a(9))e 2 ], 6 is the Heaviside step function and 
4>(v x , v y , t) is the pdf of the gas particles. The F{9) in Eq. (jSJ) is a shape function of the 
object and it is such that dl = CF(9)d9 is the length of its outer "effective" surface C 
that has a tangent between 9 and 9 + d9 (see Fig. [T]for an example). 

Also the dynamics of the pdf, <f>(v, t), of the gas particles obeys, in the dilute limit, 
a BE analogous to Eq. ([7]), with the contributions of disk-disk collisions, disk-object 
collisions and the coupling with the external driving. Then the BE for <f>(v, t) can be 
written as 

d<P(v,t) 



dt 



J[^|0,0]+ / dv' x / dv' y [W(v\v')cj>(v',t) 



W(v\v)(j)(v,t)}+B(j)(v,t). (9) 

In the above equation J[v\4>, 4>] is the Boltzmann collision operator for the disk-disk 
interactions [12] , while B is an operator representing the effects of the viscous force and 
of an external bath allowing the granular gas to reach the steady state. If we consider, 
as thermostat, the bath of Eq. ([6]), this operator has the form 

B(f)(v, t) = t)] + Q A v t)} (10) 

where A„ is the Laplacian operator with respect to the velocity 
The transition rate in Eq. ([9]) is given, instead, by 

/*2"7r ^+00 

W{v\tf) = / n ob CF(9)d9 / dVP{V',t){V' - 0) ■ k ± 



G[(V'-v>)-k\]5\v x -v' x Smf> 



e 2 K(9,a(9)) 
sin 9{V - v' x ) + v' y cos e\Y{v y ~v' y + 



cos 9 



■ [sine(V -v' x )+v' y cos9]\ (11) 



e 2 K(9,a(9)) 

where n b = 1/Ais the density of objects in the system. The function CF(9) is the same 
of Eq. (jSJ) because it is connected to the differential cross section that is a property of 
the colliding couple. 



Models of granular ratchets 



5 



A first quantitative information about this system is the collision frequency uo r c d between 
ratchet and disks. It can be obtained from the transition rate (jSJ), using the relation 



.rd 



dVP(V) / W(V'\V)dV. 



(12) 



The value of uj r c d is then determined by the choice of P(V). In our previous article [7] 
we have shown that, if the ratchet mass is comparable or smaller than the disk mass, 
i.e. e 2 > 1, the ratchet pdf is asymmetric and deviates strongly from a Maxwellian 
distribution. It is essential, in this case, to include also the third moment of the 
distribution [14J. Then we assume, for the object, a P(V) of the form 

M(V-(V}) 2 - 



P(V) 



M 
2nT r 



6dV 3 



cxp 



2T r 



(13) 



where T r = M((V — (V)) 2 ) is the granular temperature of the ratchet and £ = 
((V — (^)) 3 ) is a measure of the asymmetry of P(V) about the average value. Assuming 
that (V) <C Vth = \/2T r /M and retaining only the terms of first order in (V), we can 
perform the integrations in Eq. (fT2l obtaining that 



v. 



rd 



V^ T rd 



uAz) 



(14) 



where r T ° d 



(nCVth)' 1 , z = l/(r]e 2 ) and 77 = T r /T g is the ratio between the temper- 
ature of the tracer and the gas. Obviously £ must be small enough to make positive 
formula [TU otherwise the assumptions must be revisited. The coefficients in Eq. (I14p . 
which depend on z and on the shape of the object, are given explicitly in the Appendix. 



3. Evolution of moments 

The mean velocity and granular temperature of the object, as well as those of the 
be calculated from the two BE above, Eqs. d7j) and <^j. Starting from these, 
in fact, we obtain a set of equations for the first moments of the distributions, using 
suitable approximations to close the set. To this aim we can assume, reasonably, a 
Maxwellian distribution of the disk velocities since the gas is directly coupled to the 
bath, i.e. <p(v) = m/(2nT g ) exp[— m(v — (v)) 2 /{2T g )\. The granular temperature of the 
gas is therefore T g = m((v — (v)) 2 )/2. The equations for the first three moments of the 
distribution of P(V) can be obtained multiplying both sides of ([7j) by V, M(V — (V)) 2 
and (V — (V)) 3 respectively, and performing the integrations. 

The analogous equations for (v) and T g can be extracted in the same way, starting from 
Eq. ([9]) and considering also the contribution deriving from disk-disk collisions. 
After long calculations, with the further assumption (V) <C V t h = ^/2T r /M, we have 
derived the following equations for the moments of the object 

d(V) e 2 \V th a 2 (z) 
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a 3 (z) 



71 



3^V 2 h \ 



dT r 
~dt 



(%> + 

-^6 1 (^) + ^6 2 (^)((V r )-K))H 



M 

^h(z)(Vy) + h£ + ((V) - (V x )) + 



i/i 



at 



C3( ^-^(^-^c 5 (,)(^} + 



\/ 2 
2^ 
V 2 



V t 



th 



V 2 



C7 



W 71 v th 



V 



til 



and for those of the gas 

d(v x ) 



dt 



Nr? d L 4 



V th , x , 02(2; 
aim + 



03 U J 



7T 



(%} + 



at 



V 



th 



dn(z) 



7[ 

m 



(%> + 



I/ 2 



di(z) + 

3v^r 



«V) - (v x )) + 

71 

- T(v x ) 

((V) - (v x )) + 



th 

a 3 (z 



71 



dt 



2^ 



((VO - <*«» + 



3^ 2 



s^ 2 



2r dd 



1 — a 



th 

2 1 



2^ + 4Q 



6 

(15) 



(16) 



(17) 



(19) 



Hh 



(20) 



is the 



where v t h = \f^T g /m is the thermal velocity of the disks, r dd = (yflimavt. 
collision time between two disks and it is important to remind that z is time dependent. 
The coefficients of the above equations are given in the Appendix. Comparing the Eqs. 
(USD and flUD we obtain that 

d(v x ) 1 d(V) 



-T{v x ). 



(21) 



dt Ne 2 dt 

In the steady state this implies that (v x ) = 0. Moreover if the object is symmetric (both 
in shape and inelasticity) with respect to the axis x, then the coefficients u 3 , d\, d 3 , a%, 
b 3 , be, cq and c$ vanish (see Appendix). From Eq. ffT9]) it results that in this case the 
stationarity implies also (v y ) = 0. The average velocity of the object, (V), in the steady 



Models of granular ratchets 



7 



state is given by 

(V) = -Vjf ?iM + 3^44l. (22) 
L 4 a 2 (z) V t s h a 2 (z)l 

As previously said the contribution of £ is decisive if e 2 > 1, while it vanishes if the 
ratchet mass is very large respect to the disk mass. In this case the average velocity is 
determined by the ratio ai(z) / a 2 (z). To give an example, if we consider an exact isosceles 
triangular ratchet with angle opposite to the base equal to 26*o and with a± = a 2 = a, 
we have for e 2 < 1 

(v) _ ev th y^ j^deF{e)^e _ 

1-77 ./2^ e 2 (1 _ sin Q) (23) 



4 V m 

This is just the equation for (V) obtained in [5]. On the contrary, if a flat "piston" 
(perpendicular to the x axis) with the two faces with different inelasticities a%, for the 
left face, and a 2 , for the right face, one retrieves [7j: 



) = - \l~^ Tn~r~~ > [i + «■(■? -i)l P*) 

m 4(2 + ot dx + a sx ) L J 

which, in general, has a larger signal-noise ratio (V) / \fTg-, with respect to the uniformly 
inelastic case (1231) . and should be easier to be observed in experiments. In conclusion, 
the asymmetry of the system, appearing in the coefficients a 1; a 2 , and a 4 , determines 
a net drift of the object: a study of these coefficients shows that a modulation of the 
inelasticity along the surface, i.e. a non-constant a{9) is more efficient in producing 
a net drift, with respect to the geometrical asymmetry. It is also interesting, looking 
at Eqs. (|r5l) -( T20l) . to discuss the degree of coupling between the gas, the object and 
the thermal bath, which is determined by the three characteristic times present in the 
system: the relaxation time of the thermal bath Tb, the disk-disk collision time Tdd and 
the disk-ratchet collision time r^ r = Nr® d . Here, we are interested in the case tm < U, 
where inelastic collisions among the gas particles become relevant; in this case, three 
scenarios can occur: (i) when r^ r < Tdd < U, (ii) when Tdd < T dr < U or (iii) when 
T dd < T b < T dr- I n case (i), the dynamics of the ratchet and the disks are strongly 
coupled: in this case we expect the region of gas surrounding the object to be correlated 
with the object itself, making doubtful the assumptions of homogeneity and diluteness 
introduced at the beginning to treat the system with Molecular Chaos. If, instead, the 
disk-disk collisions are more frequent, i.e. in cases (ii) and (iii), fast and homogeneous 
relaxation of the gas is expected. In particular, in situation (iii) the gas dynamics is 
dominated by internal collisions and by the external driving, and can be regarded as 
uncoupled from the ratchet: in this limit the gas velocity pdf 4>(v) is constant and 
Eq. (CO) is a Master Equation, i.e. the ratchet velocity performs a Markov process [15J. 
The condition for the occurrence of case (i) can be approximated by 

r^w 1 I Tin fo~7 . . 
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where pd = 4m/(iTa 2 ) and p r = AttM/C 2 are proportional to the mass densities of the 
disks and the ratchet, respectively. Then case (i) occurs when pd/p r > A 2 and this 
never takes place, in practice. Cases (ii) and (iii), instead, are characterized by the ratio 
T dr/ T b being smaller or larger than 1, respectively. Because T% r oc A/C and A ^> C, 
in order to neglect the excluded area of the object, the case (ii) is obtained only if t& 
is large enough. This implies very long simulations in order to obtain average values 
statistically relevant. For this reason, we have preferred to verify the Eqs. (|15p - (|20p in 
situation (iii). In principle, and in particular in cases (i) and (ii), one should verify the 
stability of the stationary state. A study of linear stability is the objective of ongoing 
research and of a future publication. 

4. A simple case: an equilateral triangle 

The aim of this section is to compare our theoretical results with DSMC numerical 
simulations, for a particular choice of the ratcheting object. Note that the choice of 
DSMC algorithm always satisfies the Molecular Chaos. No constraints are imposed on 
the velocity pdfs: therefore this comparison is a test of the many assumptions done 
about the pdfs, to obtain Eqs. ( jT5|) - (|20T) . For the sake of simplicity, we consider as 
object an equilateral triangle with different inelasticities for the left face and for the 
two right faces (see FigJT]). In particular, we choose that the coefficient of restitution is 
«2 if < 9 < 7r, and it is ot\ if tt < 9 < 2n. The collision rules (jlj) are well defined if 
the surface is smooth. We consider then a triangle with the vertex shaped by circular 
arches of radius a v /2. In this case the outer surface is C = 3L + 7r(cr + a v ) where L is 



o 


/ 




^ o 


G 

o 










*» o 


o 




L 




\. N 

4 1 




















\ o 


o 


\ 


o 


o 


o 



Figure 1. Sketch of an inelastic equilateral triangle with different coefficient of 
restitution on left and right side [a,\ and cti respectively) and surrounding by a gas 
of inelastic particles. The dashed line correspond to the effective triangle seen by a 
colliding disk. 



the "linear" side of the triangle (see Fig. [T]). The function F(9), using the symmetry of 
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,)[0--7T 



the object, becomes 
with r e = (a + c v )/2. 

Moreover, the transition rate (jSJ) can be written as 



+ -4 
C 



(26) 



??2 



2ttT„ 



(y'-V)e(V'-V)g^V,V') + 



(V -V')Q{V -V')g + {V,Y) 



(27) 



where 



d9F{9) 

sin0>O 

m sin 2 6 1 



K 



.a 2 



sin 2 



exp < — 



2T„ 



.a 2 J 



1 2 



sin 2 



(y - V) + V 



sin6»>0 



d6F{6) X(9, a 2 , V, V) 
d6F(6) X(e,ax,V : V') 



sin6»<0 



The last line of Eq. (TJHD defines the function A(0, ot(0), V, V) 
the symmetry respect to 0, the above expressions become 



(28) 
(29) 

Using the Eq. (1261 and 



or /7r \ r f n 

g.(V, V) = -^X (-, a a , V, V'J + 4 j( d0 A(0, a 2 , V, V) 

9+(V, V) = i\ (|, ai, V, y') + | /d0 A(0, ai, 7, 7') 



(30) 
(31) 



Some sections of the transition rate surface, for different values of e 2 , are shown in Fig. 
HI together with the simulation data obtained from a DSMC with T b /r dd = 137.5 and 
add = 0.9. The Figure [2] displays a very good agreement between theory and DSMC for 
the transition rate for all values of a 2 and e 2 studied. 

In Fig. [3] we show the rescaled observable (V) /Vth as a function of the right coefficient of 
restitution a 2 for different values of e 2 . The theory results in good agreement with the 
simulation data, supporting our assumptions. The trends are analogue to those obtained 
in [7], where only the inelasticity asymmetry was considered: such a similarity indicates 
that the effects, due to the different inelasticity of the ratchet, are predominant with 
respect to the geometrical asymmetry, that however can be not neglected [5]. The data 
referring to ratio rj of the granular temperatures of the system go in the same direction 
(see Fig. HJ). 



5. Conclusion 



Within the present article we have investigated the statistical properties of a specific 
non equilibrium system, a 2D object sliding along an axis and colliding inelastically 
with a granular gas coupled to a thermal bath. If the object has an asymmetric shape 
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Figure 2. The transition rate of an inelastic equilateral triangle for different 

values of the mass ratio e 2 and the coefficient of restitution 012. The panels a, c and 
e show the trend as a function of V' for 0,2 — 0.9, add = 0.9 and e 2 = 0.01, 1 and 
10 respectively. The panels b, d and f instead show the trend as a function of V for 
«2 = 0.3 and for the same values of e 2 . The symbols correspond to the simulation 
data, while the lines are obtained from Eqs. ([?f)) and (|3T)|) - (f3"Tj) . 



or possesses a non uniform inelasticity profile, one observes a net drift. It is possible, 
employing the BE, to describe theoretically the system. We first obtain the collision 
frequency between ratchet and disks. Such a quantity determines the interplay between 
the gas and the ratchet dynamics. Secondly, we have derived a system of coupled 
equations of some relevant averages of the distribution functions of the ratchet and of 
the gas, respectively, describing the time evolution of the whole system. In the limit of 
light objects (e 2 <C 1), it is fundamental to include the third moment of the ratchet in 
order to attain a satisfactory description of the ratchet behavior. We have performed 
DSMC simulations in the case of an equilateral triangle and compared the theoretical 
predictions with the numerical results and found a fair agreement: in particular, we 
have measured the transition rate of the tracer, its mean velocity and the granular 
temperature for difference e 2 values. Whereas the theory gives an explicit representation 
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of such transition rates based on the assumption of a Maxwellian velocity distribution 
for the gas particles, our simulations give a direct estimate of the same quantity. The 
good agreement with the analytic results confirms our assumptions. 



3.0 




a 2 

Figure 3. The velocity of an equilateral triangle ratchet, rescaled with its thermal 
velocity Vth = y/T r /M, as function of the coefficient of restitution a.2 for oyd = 0.9, 
Tb/Tdd = 137.5 and for different values of the parameter e 2 : 10 (circles), 1.0 (squares), 
0.1 (diamonds) and 0.01 (triangles). The symbols correspond to the simulation data 
while the lines correspond to the solutions obtained from Eqs. (Tl5|l - (l20|) . 
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0.2 0.4 0.6 0.8 1.0 

a 2 



Figure 4. The temperature ratio r\ = T r /T g as function of the coefficient of restitution 
«2 for the same cases of Fig. |3] The symbols correspond to the simulation data while 
the lines correspond to the solutions obtained from Eqs. (p~5 |) - ([20|) . 
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The coefficients of the Eq. (fT4"l) and Eqs. (|T5l) -( fT7l) can be written as 
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Ui(z) 
u 2 {z) 
u-i{z) 
ai(z) 
a 2 {z) 
a 3 (z) 
a±(z) 
b x (z) 

h(z) 
6 4 = 

h(z) 



2tt 



d9F(9)\/z + sin 2 9 
d9F(9)- Sin49 



2- 



d9F(9)- 



2n 



d9F{6) 



z + sin 2 Of/ 2 
sin 3 6 cos 6 
z + sin 2 9fl 2 

7-^^ (1 + a(g))(z + S m a g) 
1 + e z sin 6 



2tt 



dgF(g) 



sin 2 g 



2tt 



dgF(g) 



1 + e 2 sin 2 g 
sin 9 cos g 



2- 



d9F(9) 



1 + e 2 sin 2 9 

(l + a(g))sin 4 g 



(1 + a(g))\A + sin 2 g 
(l + a(^))v / 2 + sin 2 ^ 



(1 + e 2 sin 2 0)V 2 + sin 2 

/** dgF(g) (1 + a( f )£4 2 S ^ 2 2 * V* + sin 2 g [*(1 + «(g)) + sin 2 g(a(g) - 1) - 2] 
J (1 + e^sm 9) A J 



2t 



sin 6 1 



2tt 



2tt 



dOF{9)- , . 

(1 + e z sin g)^ 

(l-a(g)e 2 sin 2 g)} 

(1 + e z sin 9) z 
(l-a(9)e 2 sin 2 g)} 
sing 



dgF(g) 



3(1 + a(9)) 2 ze A sin 2 g + [l - (2 + 3a(g))e 2 sin^ 



3(1 + a(9)) 2 ze A sin 2 9 + [l - (2 + 3a(g))e 2 sin 2 



(l-a(g)e 2 sin 2 g) 2 



(1 + e 2 sin 



2 0^2 fl\2 



2tt 



sin 2 g 



■{6^ 2 (1 - a(9)e 2 sin 2 g) 2 + 2 sin 2 g ■ 



b 6 (z) 



d9F{9) ^ + ^ gin2 ^ + gin2 ^ 5/2 

12(1 + a(g))V sin 4 g + 5(1 + e 2 sin 2 9) (l + e 2 sin 2 9-4(1 + a(9))e 2 sin 2 g) 
+2 sin 4 g(l - a(g)e 2 sin 2 9) [l - (2 + 3a(g))e 2 sin 2 g] - (1 + e 2 sin 2 g) 2 ■ 
(62 2 + 5zsin 2 g + 2sin 4 g)| 

d9F(9) j—— - . 2m5/2 {6, 2 (l - «(g)e 2 sin 2 g) 2 + * sin 2 g • 

(1 + e z sin g)^ + sin 9yi l I 

"l2(l + a(g)) 2 e 4 sin 4 9 + 5(1 + e 2 sin 2 g) (l + e 2 sin 2 g - 4(1 + a(9))t 2 sin 2 g) 
+2 sin 4 9(1 - a(9)e 2 sin 2 g) [l - (2 + 3«(g))e 2 sin 2 gl - (1 + e 2 sin 2 g) 2 • 
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(6^ 2 + 5,2sin 2 + 2sin 4 0)} 
sin0 
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{(1 + a(0))Ve 6 sin 2 9 + (1 + a(9))ze 2 1 - (1 + 2a(9)) 



Cl (z)= / d9F{9)- , . 

(1 + e^sm ft')' 3 

e 2 sin 2 0] (1 - a(0)e 2 sin 2 0) - (1 - a(9)e 2 sin 2 0) 3 } 
sin 2 



j-2-K 

c 2 (z) = / d0F(0) 



(1 + e 2 sin 2 9f{z + sin 2 0) 3 / 2 



|8^ 3 e 6 sin 2 0(1 + a(9)) 3 + Qz 2 e 2 [l + a(0) 



(1 + a(0))(l + 3a(0))e 2 sin 2 + (1 + a(0))(2 + 5a(0) + 4a 2 (0))e 4 sin 4 



+ 



j-2-K 

c 3 (z) = / d9F(9) 
Jo 



(1 + 2a(0))e 2 sin 2 9 - lj +2 sin 2 0(1 - a(0)e 2 sin 2 0) [(3 + 4a(0))e 2 sin 2 - 
1] - (1 + e 2 sin 2 0) 3 (3z + 2 sin 2 0) } 

7T— — 2 * - . 2m5/2 {2^ 3 (l - a(0)e 2 sin 2 0) 2 [l - (3 + 4a(0)) • 
(1 + e^sm 0) d (z + sm 0) t, /" : L 

e 2 sin 2 0] + z 2 sin 2 [l5 - 5(1 + lOa(0))e 2 sin 2 + 5(1 + 4«(0) - 12a 2 (0)) • 



e 4 sin 4 + (1 - 2a(9) - 12a 2 (0) - 24a 3 (0))e 6 sin 6 



+ 2z sin 4 



10- (1 + 



2^ 3 + I5z 2 sin 2 + 



c 4 



c 5 (^) 



2tt 



rf0F(0) 



31a(0))e 2 sin 2 + (1 + 4a(0) - 33a 2 (0))e 4 sin 4 - a(0)(l + 3a(0) + 12a 2 (0)) 
e 6 sin 6 0j + 8 sin 6 0(1 - a(0)e 2 sin 2 0) 3 - (1 + e 2 sin 2 0) 3 

2(ksin 4 + 8sin 6 0] | 
sin0 



(1 + e 2 sin 2 9 f 



2tt 



d9F(9) 



(1 - a(0)e 2 sin 2 0) 2 [1 - (3 + 4a(0))e 2 sin 2 0] 
e 2 sin 2 



(1 + e 2 sin 



2 o^ 2 ^3 



2M3/2 



+ sin 2 0) 



{82 3 (1 + a(0)) 3 sin 2 + 62 2 (1 + a{9)) 



c 6 (z)= / d9F(9) 
Jo 



1 - (1 + 3a(0))e 2 sin 2 + (2 + 5a(0) + 4a 2 (0))e 4 sin 4 + 6z(l + a(0)) sin 2 
3 - 6a(0)e 2 sin 2 + (1 + 2a{9) + 4a 2 (0)e 4 sin 4 0+2 sin 4 [6(1 + a{9)) + 

3(1 + «(0))(1 - 3a(0))e 2 sin 2 + (1 + 3a 2 (0) + 4a 3 (0))] } 

e 2 sin0cos0 r „. .„ N .o 9 „ 

|8^ 3 (1 + a(0)) 3 sin 2 + 6z 2 (l + a(0)) 



c 7 



2tt 



d0F(0) 



(1 + e 2 sin 2 0) 3 (z + sin 2 0) 3/2 
1 - (1 + 3a(0))e 2 sin 2 + (2 + 5a(0) + 4a 2 (0))e 4 sin 4 0] + Qz(l + a{9)) sin 2 

3 - 6a(0)e 2 sin 2 + (1 + 2a(9) + 4a 2 (0))e 4 sin 4 0+2 sin 4 [6(1 + a(9)) + 

3(1 + a(0))(l - 3a(0))e 2 sin 2 + (1 + 3a 2 (0) + 4a 3 (0))] } 
sin 



(1 + e 2 sin 2 0) 3 
« 2 (0)(3 + 4a(0))e 6 sin 6 



1 - 3(1 + 2a(0))e 2 sin 2 + 3a(0)(2 + 3«(0))e 4 sin 4 
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"l - 3(1 + 2a(0))e 2 sin 2 + 3a(0)(2 + 3a(0))e 4 sin 4 9 - 



c 8 



d 2 (z) 



2tt 



v l + e 2 sin 2 9f 
a 2 (0)(3 + 4a(0))e 6 sin 6 
cos0 



d6F{6) 



1 + e 2 sin 2 



1 + a(0)) (z + sin 2 0) 



2tt 



Jo (1 + e 2 sin 2 0)v^+ sin 2 

ei (z) = / * d9F(9) 1 + a ( 9 ] y/ z + sm 2 e \z{a{9) - 1 - 2e 2 sin 2 0) + (1 + a(0))sin 2 0l 
J Q 1 + e 2 sin 9 

/2tt 
d0F(0)- 



1 + e 2 sin 
1 + a(0))sin 3 0cos0 



sin 



2 cir, 2 fl^2 



+ e 2 sin 



3(1 + a(0)) 2 sin 2 9 + z 2(1 + e 2 sin 2 0) 2 + (1 + a(0)) 



(3a(0) - 1 - 4e 2 sin 2 0)] | 



p2tt 

e 3 (z) = j( d9F(9)- 



COS0 



e 4 = 
e 5 (^) 

ee(z) 



2tt 



rf0F(0) 



+ e 2 sin 2 0) 2 
(3a(0)-l-4e 2 sin 2 0)] } 
l + a(0)) 2 sin 3 



3(1 + a(0)) 2 sin 2 + z 2(1 + e 2 sin 2 0) 2 + (1 + a(0)) 



2tt 



d9F{9) 



(1 + e 2 sin 2 0) 2 

(1 + a(0))sin 4 



2tt 



(1 + e 2 sin 2 0) 2 (z + sin 2 0) 5 /2 
(4 + 3a(0) + e 2 sin 2 9) + 3(1 + a{9)) sin 4 
(1 + a(0))sin 3 0cos0 



z\h + 3a (0) + 2e 2 sin 2 0) + 2z sin 2 



d0F(0)- 



(1 + e 2 sin 2 0) 2 (z + sin 2 0) 5 / 2 
(4 + 3a(0) + e 2 sin 2 0) + 3(1 + a{9)) sin 4 



z 2 (5 + 3a(0) + 2e 2 sin 2 0) + 2z sin 2 
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